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An  Objective  Method  for  Forecasting  Solar  Flares 


I.  INTRODl  C’TION 

This  report  is  a  continuation  of  an  earlier  study  (Hirman  et  al,  1980)  in  which 
multivariate  discriminant  analysis  (MVDA)  is  used  in  a  computer  program  to  pro¬ 
duce  an  objective  daily  solar  flare  forecast.  The  essential  feature  of  the  statistics 
package  is  the  comparison  between  a  number  of  input  parameters  and  a  number  of 
output  classes,  in  which  the  discrimination  between  the  classes  in  terms  of  the 
input  parameters  is  maximized  by  constructing  appropriate  classification  functions. 
In  the  application  to  flare  prediction,  the  input  parameters  are  daily  solar  param¬ 
eters  for  each  active  region  on  the  solar  disk,  and  the  output  classes  are  the  levels 
of  flare  activity  occurring  the  following  day  within  the  same  active  regions.  We 
have  used  more  than  two  years  of  data,  of  which  approximately  25  percent  has 
been  used  to  derive  the  classification  functions.  The  latter  are  then  extrapolated 
forward  in  time  to  produce  a  true  forecast. 

The  computer  program,  known  as  UMD07M,  was  originally  written  at  I'CLA,  1 
although  the  particular  version  used  here  was  developed  further  by  Seagraves“  to 

Keceived  for  publication  3  Feb  1981 

1.  Dixon,  W.  J.  (ed.)  (1968)  Biomedical  Computer  Programs,  University  of 

California  Publications  in  Automatic  Computations,  No.  2,  University  of 
California  Press,  p.  214a. 

2.  Seagraves,  P.  H.  (1972)  UBC  BMD07M  Stepwise  Discriminant  Analysis, 

1  niversity  of  British  Columbia  Computing  Centre  1  locumentation. 


include  tho  Cooley  ami  l.ohnes  classification  procedure,  ’  and  the  I  .nchenbrueh 
\-l  technique.  *  The  Cooley  and  l.ohnos  procedure  does  not  assume  uniformit\ 
of  variance,  and  this  sometimes  results  in  better  classification  scores.  The 
computational  burden,  however,  is  increased  becau.-e  linear  classification 
functions  are  not  possible:  instead,  canonical  variables,  const  meted  from  the 
original  input  parameters,  are  used  as  a  transformation  t  ■  reduce  the  matrix 
dimension  in  the  classification  formulas.  1‘he  1  .achenbruch  technbjue  retr.oves 
bias  when  the  procram  classifies  its  own  data  base. 

A  complete  description  of  the  mathematics  is  beyond  the  scope  of  this  report. 
Tho  reader  may  consult  \nderson  and  Kao1’  for  references  on  discriminant 
analysis.  A  discussion  of  the  suitability  of  apply  inn  various  statistical  methods  to 
discrete  input  variables  is  contained  in  Vccchia  et  al.  '  The  latter  point  is  of 
particular  interest  because  the  work  of  Vccchia  et  al  uses  tho  same  discrete  data 
base  as  used  herein,  to  produce  solar  flare  probability  forecasts  using  discrimin¬ 
ant  analysis  (without  the  Cooley  and  l.ohnos  procedure)  and  logistic  regression 
analysis. 

An  important  feature  of  the  present  study  is  the  comparison  of  the  objective, 
computer  forecast  with  a  subjective,  conventional  forecast  prepared  during  the 
same  test  period  for  the  same  active  regions  on  the  sun.  Without  such  a  bench¬ 
mark  for  relative  evaluation,  tho  presentation  of  anv  forecast  method  has  con¬ 
siderably  reduced  merit. 

2.  1)  VT  \ 

The  data  used  herein  were  obtained  from  the  region  analysis  program  at  the 
\0.\A  Space  Environment  Services  Center  (STCSO)  in  Boulder,  Colorado.  The 
region  analysis  program  collects  daily  a  variety  of  solar  parameters  for  each 
active  region  on  the  solar  disk,  it  is  important  to  note  that  there  is  no  attempt 


8.  Cooley,  W.  W. ,  and  l.ohnes,  P.  R.  (1962)  Multivariate  Procedures  for  the 
Behavioral  Sciences,  Wiley,  New  York. 

I.  Lachenbruch,  P.  A,,  and  Mickey,  M.  R.  (1968)  Technometrics,  10:1. 

i.  Anderson,  T.  U  .  (1958)  An  Introduction  to  .Multivariate  Statistical  Analysis, 
Wiley,  New  York. 

6.  Rao,  C,  R,  (1974)  Advanced  Statistical  Methods  in  Biometric  Research, 

Hafner. 

7.  Vecchia,  D.  Caldwell,  G.  A.,  Tryon,  P.  V.,  and  Jones,  R.  H.  (1980) 

in  Sol. -Terres.  Pred.  Proc. ,  Vol.  3,  R.  F.  Donnelly  (ed. ),  C-76. 
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in  this  program  to  select  the  more  flare-productive  regions.  The  parameters 
include  radio  and  \-ray  data,  but  most  are  derived  from  optical  data  supplied 
by  the  I  SA  1'/  A \VS  SOON  system.  The  parameters  contain  information  which  the 
SESC  forecasters  consider  vital  to  the  preparation  of  a  24 -hour  flare  forecast. 
The  present  study  uses  data  for  the  period  1  January  1977  to  .11  January  1979, 
containing  fiOlVi  active-region  days  (records)  that  have  been  checked  for  errors 
and  internal  consistency.  Random  scrutiny,  however,  lias  shown  that  errors 
still  remain.  After  several  reassignments  of  parameter  values  and  definitions, 
we  arrived  at  the  form  of  the  data  base  shown  in  Table  1, 


Table  1.  SESC  Region  Analysis  Parameters  (Modified) 


PARAMETER 


ASSIGNED  V  A I  .I  E 


1.  DATE 

2.  REGION  NUMBER 

3.  REGION'S  PIRST  APPEARANCE  LONGITUDE 

4.  CURRENT  LONGITUDE 

5.  N/S  LATITUDE 

6.  CURRENT  LATITUDE 

7.  CARRINGTON  LONGITUDE 

8.  REGION  AGE 

9.  SPOT  CLASS  1 


Table  1.  SFSC  Region  Analysis  Paramoters  (Modified)  (continued) 


10.  SPOT  Cl. ASS  2 


r/x .  1 


a .  3 

h .  4 

k  . . r> 


11.  SPOT  CLASS  3 


x .  1 

o .  2 

i  .  3 

c  . . 4 

12.  MAGNETIC  CLASS 

No  spots . 0 

Alpha .  1 

Beta . . .  2 

Beta-Gamma .  3 

Gamma .  4 

Beta-Delta .  5 

Beta-Gamma-Delta .  6 

Gamma-Delta .  7 

13.  MAGNETIC  POLARITY  OF  STRONGEST  FIELD .  (+/-) 

14.  MAGNETIC  FIELD  STRENGTH . . .  (Gauss) 

15.  MAGNETIC  GRADIENTS .  (Gamma/km) 

16.  INTERACTION  WITH  ANOTHER  REGION 

None  . . . . . .  0 

Spots  of  opposite  polarity  converge  (from  less 
than  two  degrees  apart) . . .  1 

17.  SUNSPOT  DYNAMICS 

No  spots  or  no  motion .  0 

Coalescing  of  spots  . .  1 


Table  1.  SE.SC  Region  Analysis  Parameters  (Modified)  (continued) 


Spot  rotation . 2 

Relative  motion  between  oppositely  poled  spots . 

18.  STAGE  OK  DEVELOPMENT 

No  spots . 0 

Mature  group  (stable)  . . .  .  .  .  1 

Decaying . 2 

Growing . 8 

Rapid  decay  (spot  numbers/areas  decrease  by  ■»  dO  ) .  4 

Rapid  growth  (spot  numbers/ areas  increase  bv  >  dO'" ) .  f> 

Rapid  growth  (>  100"') .  6 

19.  LEADER/TRAILER  FIELDS 

Structure  not  definite . 0 

Returning  Region  . .  1 

<5  deg  of  neutral  line  and  out  of  phase . .  2 

> d  deg  of  neutral  line  and  in  leader  fields .  3 


>5  deg  of  neutral  line  and  in  trailer  fields 
< d  deg  of  neutral  line  and  in-phase . 

20.  RETt  RNING  REGION 

21.  SECTOR  BOUNDARY  RELATIONSHIP 

22.  ASSOCIATED  FILAMENT 

None . . . 

Filament  unchanged  . . 

Filament  growing . . 

Filament  disappeared  within  past  24  hrs  . 
Filament  darkens  or  is  active . 

23.  EMBEDDED  FILAMENT 

None . . . 

Filament  present . 


Active  filament 


Table  1.  SESC'  Hogion  Analysis  Parameters  (Modified)  (continued) 


Table  1.  SESC  Region  Analysis  Parameters  (Modified)  (continued) 


31.  ISOLATED  POLE 

32.  EMERGING  FLEX 

None,  or  region  is  new . . 

New'  flux  emerges  within  spot  group . 

New  flux  emerges  near  region  (within  5  deg) . 

33.  ARCH  FILAMENT  SYSTEM 

34.  RADIO  UURST/SWEEP 

Xone  occurred  . . 

>250  flux  units  at  10  cm . . 

>1000  flux  units  at  10  cm . 

Type  III . 

Type  IY . 

Type  II  and  IV  .  . . . 

F  Hurst . . . 

Major/complex  10  cm  burst . 

>1000  flux  units  at  10  cm  plus  a  U  burst,  or 
Type  III  and  IV,  or 

250  flux  units  at  10  cm  plus  Type  III  and  IV  . . . 

35.  REGION'S  FIRST  APPEARANCE  (TRANSIT  HISTORY) 

36.  FLARE  HISTORY 

No  flares  have  occurred . 

C  class  flares  have  occurred  . . . 

M  class  flares  have  occurred . 

X  class  flares  have  occurred . 

37.  FLARES  TODAY 


0 

1 


0 

1 

2 

3 

4 

5 

6 
7 


8 


0 

1 

2 

3 


None 


38.  PHOTON  HISTCin 


None  occurred . 0  j 

Proton  event  occurred . 1  j 

Ground  level  event . . . . . 

.19.  PHOTONS  TODAY  ] 

None . 0 

Occurred  .  1 

40.  REGION  F<  'HKCASTS  (SESC)  j 

l 

Probabilities  Tor  each  class  of  flare  (none,  G,  \l,  or  \)  for  1 

each  region.  Tor  the  34-hour  pnrioii  beginning  at  0  hr  1  T  next  1 

day.  Proton  event  probabilities  are  similarly  stated. 


Most  of  the  parameters  in  Table  1  have  oeen  assignee  discrete  values 
according  to  categories  which  are  subiectivelv  related  to  increasing  fl are 
activity.  This  subjectivity  is  the  weakest  link  in  an\  scheme  utilizing  oiiiective 
procedures  for  prouucing  a  forecast  solelv  from  data,  in  essence,  the  situa¬ 
tion  merely  allows  the  element  of  subjectivity  to  reside  entirely  in  the  data 
acquisition  process.  Probably,  this  situation  is  preferable  to  having  subiecti- 
vity  introduced  also  in  the  forecast  preparation.  'There  are  several  parameters 
(e.g.  spot  class,  flare  history,  magnetic  class)  for  which  assigned  values  are 
based  upon  quantitative  studies.  Fortunately,  for  perhaps  therefore;  )  these  pararr, 
eters  are  among  those  from  which  the  obiective  forecast  derives  most  of  its  skill. 

Perhaps  the  most  unfortunate  circumstance  is  that  for  a  large  number  of 
records  one  or  more  parameters  is  missing.  In  the  computer  program,  mis¬ 
sing  data  codes  are  replaced  by  averages  for  the  particular  parameter  in  the 
set  of  records  used  in  deriving  the  classification  functions.  Missing  data,  in 
addition  to  errors,  makes  the  testing  of  objective  techniques  difficult,  espec¬ 
ially  for  determining  the  relative  significance  of  various  parameters.  In 
order  to  portray  some  feeling  for  the  degree  of  representation  in  the  data  base 
we  note  the  following:  for  three  commonly  observed  parameters.  Spot  Class 
2,  'Magnetic  Class,  and  "Flares  Today,  "  only  5893  of  the  total  G09.>  records 
contain  all  three;  if  "Bright  Points,  "  "Spot  Class  3,  "  "Spot  Class  1,  "  "Magnetic 
Gradients,  "  and  "Sunspot  Dynamics"  are  added  to  the  first  three,  only  3732 
records  remain;  and  for  a  total  of  15  of  the  31  usable  parameters,  only  510 
records  contain  all  15.  This  is,  indeed,  a  hardship  for  statistical  analysis. 


Nevertheless,  we  are  able  to  show  later  that  at  least  some  of  these  frequently 
missing  parameters  contain  valuable  predictive  information. 

The  data  base  contains  daily,  region-by-region  entries  for  the  actual  flare 
activity,  in  addition  to  the  official  SESC  subjectively  derived  flare  forecast. 
Thus,  the  information  required  for  objective  forecast  testing,  as  well  as  for 
comparison  with  the  SESC  forecast,  is  contained  in  the  same  base.  Flares  are 
listed  according  to  their  peak  soft  (1-11  1)  X-ray  flux  at  1  Al  : 


Class  C: 

10- 

- 

R  <  10‘;’ 

-a 

watt 

_*> 

m 

c  lass  M ! 

10*-’ 

-4 

<; 

o 

V 

watt 

m 

_  •> 

Class  X : 

10 

< 

E 

watt 

m  “  . 

From  the  standpoint  of  geophysical  environment  studies,  the  classes  ft!  and  X 
are  of  greatest  importance. 

In  addition  to  Table  1,  six  combination  parameters  (Table  2),  derived  from 
certain  original  parameters,  were  included  as  input  parameters.  These  six  were 


Table  2.  Combination  Parameters  (Numbers  in  right-hand 
column  refer  to  original  parameter  number  in  Table  1) 


i 

New  Parameter  No. 

Parameter  Formula 

1 

9*  10*  11 

2 

9.  10*  IP  12 

3 

9*  10.  11*32 

4 

14.  15*  (17+25) 

5 

12-  17.  27 

6 

17-  (25+27+28) 

found  to  have  possible  predictive  significance  in  the  earlier  study  where  twenty 

Q 

such  combination  parameters  were  tested.  The  derivation  of  combination  param¬ 
eters  is  based  on  intuitions  about  the  form  in  which  predictive  information  might 
be  contained  in  the  data,  and  about  physical  quantities  (e.  g.,  energy  stored  in 
sheared  magnetic  fields)  presumed  relatable  to  flares.  The  subject  of  these  and 
other  combination  parameters  will  be  discussed  in  a  later  section. 


8.  Hirman,  J.  W.,  Neidig,  D.  F.,  Seagraves,  P.  H.,  Flowers,  W.  E.,  and 

Wiborg,  P.  H.  (1980)  in  Sol. -Terrest,  Fred.  Proc.,  Vol.  3,  H.F.  Donnelty 
(ed.),  C-64. 
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Tilt'  region  analysis  parameters  for  today  are  independent  of  any  information 
on  flare  activity  occurring  tomorrow;  therefore,  they  can  lie  used  in  practice, 
today,  to  produce  a  Hare  forecast  for  tomorrow,  assuming  that  predictive  informa¬ 
tion  is  present  in  the  parameters.  We  have  used  the  first  \  records  (with  X  =  1500, 
as  described  below)  as  a  training  set'  in  order  to  derive  the  classification  func¬ 
tions  for  three  possible  outcomes:  "No  flare,  "  V  flare,  "  and  "M  or  X  flare.  " 

M  and  \  Hares  were  grouped  together  as  a  single  class  in  order  to  reduce  statisti¬ 
cal  noise  caused  by  the  relatively  few  cases  of  larger  Hares.  The  classification 
functions  were  then  applied  to  new  records,  using  only  the  input  parameters,  in 
order  to  produce  a  true  forecast.  The  latter  procedure  was  accomplished  in 
steps  of  250  records  each,  with  the  training  set  sliding  forward  in  time,  2,50 
records  (approximately  one  month)  after  each  step.  Thus,  for  a  150f)-rocord 
training  set,  the  remaining  •!  >9 5 -record  test  set  require  If!  individual  subtests 
of  250  records  each  (except  for  the  nineteenth).  This  sliding  base  technique  main¬ 
tains  a  constant  X  records  in  the  training  set,  thereby  assuring  that  the  program 
is  trained  on  recent  data  relative  to  the  test  subset.  This,  combined  with  the 
relatively  small  size  of  the  test  subset,  minimizes  the  effects  of  secular  trends, 
either  of  observational  or  solar  origin,  which  might  be  present  in  the  data. 

The  computer  program  was  trained  on  the  X-ray  class  of  the  largest  event 
(No  flare,  (.'  blare,  or  M  <v  X  flare)  occurring  in  the  region  in  the  24 -hour 
period  following  the  acquisition  date  of  the  input  parameters.  Thus,  the  computer 
forecast  is  expressed  in  terms  of  probabilities  for  the  largest  event  to  be  in  one 
of  these  classes.  The  outcomes  are  mutually  exclusive,  with  the  sum  of  probabili¬ 
ties  over  all  classes  equal  to  unity.  The  s’ ESC  forecast,  however,  is  a  probability 
forecast  for  the  occurrence  of  each  class  of  event;  i.e.,  a  non-exclusive  format. 

In  order  to  assess  the  quality  of  the  computer  forecast,  we  derived  a  comparison 
forecast  in  the  exclusive"  format  by  selecting  the  largest  event  class  in  the 
SESC  forecast  that  was  assigned  a  probability  greater  than  or  equal  to  0.  5.  Al¬ 
though  this  is  not  an  SESG  forecast,  it  is  probably  representative  of  what  would 
be  extant  if  the  SESC  chose  to  cast  their  predictions  in  this  mode. 

In  the  following  test  results  we  present  the  forecasts  according  to  both  the 
standard  multivariate  discriminant  analysis  (MVDA)  and  the  Cooley  and  I.ohnes 
procedure  (MVDA/CL),  There  are  important  differences  in  the  character  of 
these  two  forecasts,  which,  as  will  be  shown  later,  may  be  used  to  advantage. 
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1.  TEST  KESl  I  IS 


1.1  hi' 1 1 ill m .1 1  \  Discussion 

We  aro  concerned  mainly  with  the  behavior  of  the  computer  forecast*  rela¬ 
tive  to  the  comparison  forecast  when,  for  example,  changes  are  made  m  the 
size  of  the  training  set,  choice  of  input  parameters,  solar  activity  levels,  and 
percent  of  missing  data.  In  all  cases  we  present  the  computer  forecast  along 
with  the  comparison  forecast  for  the  same  set  of  test  records.  Also  included  is 
a  list  of  input  parameters  submitted  to  analysis,  along  with  their  frequency  of 
selection  in  classifying  the  three  outcomes,  .Note,  however,  that  due  to  the  250- 
reeord  increment  the  training  sets  are  independent  of  each  other  only  when  sepa¬ 
rated  by  .-tx  or  more  subsets. 

\s  a  first  step,  we  eliminated  11  parameters  which  were  not  selected  in  any 
of  the  Id  subsets,  following  this,  the  program  was  run  again  using  the  remain¬ 
ing  2d  input  parameters.  The  results  are  given  in  Tables  A,  4,  and  a.  This 
test  ( \)  will  ,-erve  as  an  example  for  the  displays  used  elsewhere  in  tills  report. 

Table  :t  show.-  the  actual  matrix  of  region -day  forecasts  vs  r<  gam-day 
largest  events,  for  the  three  forecasts.  Table  5,  derived  from  the  data  in  Table 
1,  summarizes  the  following: 

T  Percent  of  forecasts  correct  in  the  given  event  class 
E  Percent  of  region-day  largest  events  which  were  forecasted 

\  (P  s  K)/j 

P  Plimatology  (percent  of  the  total  number  of  events  in  the  class) 
l  l  nweighted  mean  of  the  A's  for  all  three-event  classes 
W  Weighted  mean  forecast  accuracy  (the  sum  of  the  matrix 

diagonal  elements  divided  by  the  total  number  of  forecasts, 
or  events,  ill  all  classes) 

Off  1  Percent  of  forecasts  that  are  one  matrix  element  away  from 
the  diagonal 

off  2  Percent  of  forecasts  that  are  two  matrix  elements  away  from 
the  diagonal 

’These  various  scores  are  of  interest  because  of  the  several  ways  in  which 
forecasts  can  he  used.  For  example,  the  F  sco rp,  or  percentage  of  forecasts 
that  are  correct,  is  the  quantity  of  interest  to  a  customer  who  cannot  tolerate 
false  alarms.  A  quite  different  requirement  applies,  however,  in  a  situation 
where  surprise  flares  are  unwelcome.  In  the  latter  case,  the  E  score  is  the 
important  measurement.  Of  course,  knowing  the  customer's  need  in  advance 
allows  the  forecast  to  be  biased  either  toward  underprediction,  which  tends  to 
improve  the  1*  score,  or  toward  overprediction,  which  improves  the  K  score. 
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As  a  measure  of  the  "balanced"  accuracy  of  a  forecast  in  a  given  event  class 
we,  therefore,  introduce  the  average  of  K  and  E,  given  by  A. 

The  accuracy  of  a  forecast  is  always  dependent  upon  the  climatology  for 
the  event  being  forecasted.  Higher  climatological  probabilities  tend  to  improve 
the  chances  for  predictions  to  be  correct.  For  example,  it  is  easy  to  predict 
"No  Flare"  with  90  percent  accuracy,  simply  because  no  flare  occurs  in  almost 
90  percent  of  all  active-region  days.  In  comparing  cumulative  scores  between 
forecasts  it  is  imperative  to  note  the  climatology  which  prevailed  during  the 
test  period.  Climatology  is  affected  bv  a  number  of  factors,  including  event 
classification  criteria,  duration  of  forecast  interval,  and  level  of  solar  activity. 
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Table  4.  Parameters  Submitted  to  Analysis  and 
Their  Frequency  of  Selection  in  19  Subsets--Test  A 


Flares  Today 

19 

New  No.  1 

9 

Mag.  Pol. 

5 

Bright  Points 

19 

Mag.  Grad. 

9 

Neut.  I..  Chg. 

5 

New  No.  2 

17 

Mag.  Class 

G 

Spot  Class  3 

4  ! 

Spot  Dvnam. 

12 

Radio  B/S 

G 

Spot  Class  2 

3  i 

New  No.  a 

12 

Flare  Hist. 

G 

Spot  Inter. 

3 

Proton  Hist. 

1  1 

New  No.  .1 

G 

Emerg.  Flux 

1 

i  ! 

Spot  Class  1 

9 

New  No.  4 

6 

j 

Table  5.  Comparison  of  Forecast  Scores-  - Test  A 


[•‘or  pc  aster 

Event 

F 

E 

A 

C 

I 

\\ 

(  O'f  1 

(  Iff  2 

COMPARISON 

No  Flare 

93.4 

85.4 

89.4 

86.  1 

_ 

C 

2G.  2 

40.3 

33.  3 

10.  8 

52.8 

79.  2 

18.  G 

2.  2 

MAX 

28.  1 

42.  9 

35.  6 

3.2 

MVD  A 

No  Flare 

9.3.  9 

34.7 

89.  3 

86.  1 

(’ 

26.  8 

41.7 

34.  3 

10.8 

53.  6 

78.  9 

18.  4 

2.  6 

MX  X 

26.  9 

47.6 

37.3 

3.  2 

MAD  A/Cl. 

No  Flare 

91.  1 

94.6 

92.8 

86.  1 

C 

37.7 

28.7 

33.2 

10.8 

54.  3 

85.  5 

12.8 

1.7 

MAX 

41.6 

32.  0 

36.  8 

3.2 

17 


In  essence,  climatology  is  directly  dependent  upon  "bin  size,  "  Failure  to  state 

climatological  conditions  clearly  (an  unfortunately  common  practice)  makes 

9 

intercomparison  of  forecasts  almost  impossible.  It  seems  that  this  point  can¬ 
not  be  emphasized  enough. 

Because  "No  Flare"  constitutes  the  majority  of  situations  on  the  sun,  it  comes 
as  no  surprise  that  solar  flare  forecasts  are  usually  quite  accurate  overall;  i.e. , 
their  weighted  means  (\V)  are  high.  It  is  of  greater  interest,  however,  to  predict 
flares  than  quiet  conditions  and,  for  this  reason,  the  unweighted  score  1,  given 
simply  by  the  mean  of  the  A  scores  over  all  classes,  has  been  included  in  Table  a. 

Finally,  we  note  that  if  a  forecast  is  in  error,  it  is  better  to  be  wrong  by  one 
event  class  than  by  two.  Thus,  the  tendency  for  the  off-diagonal  entries  in  the 
matrix  to  cluster  near  the  diagonal  is  an  important  measure  when  comparing  fore¬ 
cast  scores  which  are  similar  otherwise.  Table  5  includes  a  measure  of  this 
error  distribution  in  the  form  of  the  Off  1  and  Off  2  scores. 

The  scores  (F,  E,  and  A)  have  uncertainties  of  approximately  ±1,  ± 3 ,  and 
±5  for  No  Flare,  C  Flare,  and  M  &  X  Flare,  respectively.  The  F  and  W  scores 
have  uncertainties  of  about  ±1.  Thus,  in  terms  of  A  and  F,  the  three  forecasts 
in  Table  5  are  essentially  identical.  The  MVDA/CL  forecast  definitely  excels  in 
the  VV  score,  although  this  is  mainly  due  to  its  tendency  for  underprediction,  which 
places  a  large  number  of  forecasts  in  the  No  Flare  column.  The  tendency  for 
underprediction  in  the  MVDA/CL  forecast  is  evident  also  in  the  F  scores  for  C, 
and  M  &  X  flares,  being  significantly  higher  than  the  corresponding  E  scores. 

On  the  other  hand,  both  the  comparison  and  the  MVDA  forecast  are  biased  toward 
overprediction.  Their  overall  similarity  is  quite  striking. 

4.2  Effect  of  Training  Set  Size 

The  number  of  records  to  be  used  in  the  training  set  should  be  large  enough 
to  provide  sufficient  statistics  to  train  the  computer  program,  yet  small  enough 
to  avoid  the  effects  of  trends  in  the  data.  The  optimum  number,  while  not  known 
from  theory,  may  be  determined  empirically  by  varying  the  training  set  size  and 
comparing  the  scores  of  the  resulting  forecasts.  Table  6  shows  the  results  for 
training  sets  of  750  and  2095  records.  Together  with  Table  5  (1500-record  train¬ 
ing  set)  we  find  differences  of  only  small  significance.  A  close  examination  of 


9.  Simon,  P. ,  Smith,  J.  B.,  Ding,  Y. ,  Flowers,  W. ,  Guo,  Q. ,  Harvey, 

K.  L. ,  Hedeman,  R.,  Martin,  S.  F. ,  McKenna  Lawlor,  S. ,  Lin,  V., 
Neidig,  D. ,  Obridko,  V.  N. ,  Dodson  Prince,  H.,  Rust,  D. ,  Speich,  D. , 
Starr,  A.,  and  Stepanyan,  N.  N.  (1980)  in  Sol. -Terres.  Pred.  Proc., 
Vol.  2,  R.  F.  Donnelly  (ed.),  p.  287. 
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Table  6.  Comparison  of  Scores  Using  750-Record 
and  2095-Kecord  Training  Sets--Test  1$  and  (' 


Forecaster 

Event 

F 

E 

A 

U 

VV 

Off  1 

Off  2 

MVDA 

No  Flare 

93.9 

86.0 

89.  9 

750  Records 

C 

28.  4 

41. 3 

34.  8 

53.  3 

80.0 

17,0 

3.0 

M  &  X 

25.  2 

45.  1 

35.  1 

MVDA/CL 

No  Flare 

91.  7 

93.  7 

92.7 

750  Records 

C 

36.  2 

32.7 

34.4 

52.  8 

85.  1 

13.  1 

1.8 

M  &  X 

36.3 

26.0 

31.  2 

MVDA 

No  Flare 

93.  9 

84.0 

89.0 

2095  Records 

C 

25.  9 

42.8 

34.4 

53.  6 

78.  4 

19.4 

2.  2 

M  &  X 

28.  6 

46.4 

37.  5 

MVDA/CI. 

No  Flare 

91.0 

95.0 

93.0 

2095  Records 

C 

38.0 

29.4 

33.  7 

54.  3 

85.  8 

12.8 

1.4 

M  &  X 

45.8 

26.4 

36.  1 

the  trend  in  the  various  scores,  however,  suggests  that  there  may  be  some 
improvement,  especially  in  the  MVDA/CL  forecast,  as  the  size  of  the  training 
set  is  increased  from  750  to  1500  records.  The  improvement  is  less  certain  in 
increasing  the  set  from  1500  to  2095.  According  to  motivations  which  will  be 
described  later,  the  E  score  is  of  interest  in  the  case  of  the  MVDA  forecast, 
while  the  F  score  is  of  prime  importance  for  MVDA/CL,  Noting  these,  the  U 
scores,  and  the  fact  that  we  do  not  wish  to  make  the  training  set  unnecessarily 
large,  we  have  decided  to  use  1500  records  in  all  training  sets. 

4.3  Inclusion  of  Additional  Combination  Parameters 

Table  4  indicates  that  five  of  the  six  combination  parameters  from  Table  2 
were  retained  for  analysis  after  the  initial  parameter  selection.  Because  several 
of  these  ranked  highly  in  frequency  of  selection  in  Test  A,  we  decided  to  test 
additional  combination  parameters.  As  in  the  case  of  the  original  six,  the  addi¬ 
tional  parameters  were  derived  on  the  basis  of  intuition.  Their  formulas  are 
given  in  Table  7. 

The  20  new  combination  parameters,  in  addition  to  the  20  parameters  used 
in  Test  A,  were  submitted  to  analysis  in  Test  D  (Tables  8  and  9),  It  is  convenient 
to  defer  the  discussion  of  the  latter  to  the  following  section. 


Table  7.  Additional  Combination  Parameters  (Numbers 
in  right-hand  column  refer  to  original  parameter  numbers 
in  Table  1) 


New  Parameter  No. 

Parameter  Formula 

Rates  of  Change 

7 

29  (today)  -  29  (yesterday) 

8 

37  -  37 

9 

(9-  10-  IP  12)  -  (9-  10-  IP  12) 

10 

17  -  17 

11 

(12-  17-  27)  -  (12*  17-  27) 

12 

38  -  38 

13 

9-9 

14 

(9-  10-11)  -  (9-  10-11) 

15 

15-15 

16 

12  -  12 

Parameters  Squared 

17 

292 

18 

372 

19 

(9-  10-  IP  1 2) 2 

20 

172 

21 

92 

22 

(New  7)2 

23 

(New  8)2 

24 

(New  9>2 

25 

(New  10)2 

26 

(New  13)2 

5>n 


Table  8.  Parameters  Submitted  to  Analysis  and  Their  Frequency 
of  Selection  in  19  Subsets--Tests  D,  E,  F,  G,  and  H 


No.  of 

Test 

Parameters 

D 

40 

Flares  Today 

19 

Hadio  B/S 

6 

Flare  Hist. 

2 

New  18 

17 

New  9 

6 

New  4 

2 

New  2 

16 

New  12 

6 

New  23 

2 

Bright  Pts. 

14 

New  1 

5 

Spot  Class  2 

1 

New  19 

12 

Neut.  L.  Chg. 

5 

Emerg.  Flux 

1 

New  1 5 

10 

New  5 

5 

New  17 

1 

Mag.  Grad. 

9 

New  14 

5 

New  19 

1 

Proton  Hist. 

9 

New  21 

5 

Spot  Class  1 

0 

New  3 

9 

New  22 

5 

New  11 

0 

New  8 

9 

Mag.  Pol. 

4 

New  13 

0 

New  20 

8 

New  7 

4 

New  16 

0 

Mag.  Class 

7 

Spot  Inter. 

3 

New  26 

0 

New  10 

7 

New  25 

3 

Spot  Class  3 

6 

New  20 

2 

a 

20 

See  Table  4 

j 

E 

15 

Flares  Today 

19 

Spot  Class  3 

14 

Hadio  B/S 

6 

Bright  Pts. 

19 

Spot  Dynam. 

13 

Spot  Class  1 

5 

Mag.  Class 

16 

Proton  Hist. 

11 

Mag.  Pol. 

Mag.  Grad. 

16 

Flare  Hist. 

10 

Emerg.  Flux 

4 

Spot  Class  2 

14 

Neut.  L.  Chg. 

6 

Spot  Inter. 

3 

F 

8 

Flares  Today 

19 

Spot  Class  2 

16 

Spot  Dynam. 

11 

Bright  Pts, 

19 

Spot  Class  3 

14 

Spot  Class  1 

5 

Mag.  Class 

17 

Mag.  Grad. 

13 

G 

5 

Flares  Today 

19 

Mag.  Class 

18 

Spot  Class  3 

15 

Bright  Pts. 

19 

Spot  Class  2 

16 

H 

3 

Flares  Today 

19 

Mag.  Class 

19 

Spot  Class  2 

19 

21 


Table  9.  Effects  of  Reduction  in  the  Number  of  Input  Parameters 


Forecaster 

Number  of 

Parameters 

U 

W 

Off  1 

Off  2 

R 

COMPARISON 

52.8 

79.2 

18.  6 

2,  2 

2.  22 

TEST  D 

MVDA 

40 

53.8 

79.  5 

18.4 

2.  1 

2.  38 

MVDA/ CL 

54.  6 

85.  1 

13.4 

1.  5 

0.73 

TEST  A 

MVDA 

20 

53.  6 

78.  9 

18.4 

2.  6 

2.  63 

MVDA/ CL 

54.3 

85.  5 

12.  8 

1.7 

0.60 

TEST  E 

MVDA 

15 

52.4 

78.  1 

18.  7 

3.  2 

2.80 

MVDA/ CL 

53.9 

85.4 

12.8 

1.  8 

0.65 

TEST  F 

MVDA 

8 

52.2 

78.0 

18.6 

3.4 

2.83 

MVDA/ CL 

53.  3 

85.0 

13.4 

1.6 

0.71 

TEST  G 

MVDA 

5 

53.0 

77.  G 

18.  5 

3.  9 

3.02 

MVDA/ CL 

53.7 

84.4 

14.0 

1.7 

0.83 

TEST  H 

MVDA 

3 

51.  5 

78.2 

17.  5 

4.3 

2.  72 

MVDA/ CL 

53.  5 

85.2 

13.  3 

1.  5 

0.  61 

•1.4  Reduction  in  the  IN  umber  of  Parameters 

The  computer  forecast  was  subjected  to  a  series  of  reductions  (Tests 

E,  F, 

G,  and  H)  in  the  number  of  input  parameters,  according  to  Table  8,  with  the 
corresponding  forecast  results  summarized  in  Table  9.  Table  9  displays  the 
effects  of  parameter  reduction  beginning  with  40  parameters  and  ending  with 
only  three.  In  addition  to  the  previously  used  scores  we  introduce  R,  the  ratio 
of  the  number  of  matrix  entries  below  the  diagonal  to  the  number  above  the 
diagonal.  This  ratio  provides  a  measure  of  the  asymmetry  of  the  forecast,  with 
values  greater  than  unity  indicating  overprediction,  and  values  less  than  unity 
indicating  underprediction. 

Table  9  clearly  illustrates  that  the  reduction  in  the  number  of  parameters 
has  a  small  but  unfavorable  effect  on  the  computer  forecasts.  We  may  regard 
the  tendencies  for  R  to  depart  further  from  unity,  for  Off  2  to  increase,  and  for 
U  to  decline,  as  evidence  for  progressively  worsening  forecasts.  These  three 
effects  are  most  noticeable  in  the  MVDA  forecast,  while  the  latter  effect  alone  is 
marginally  evident  in  MVDA/ CL. 


The  effects  of  the  parameter  reduction  are  offset  by  the  increase  in  the 
number  of  records  containing  all  or  most  of  the  parameters  submitted  for  anal¬ 
ysis  in  the  reduced  sets.  This  improvement  in  representation  occurs  because 
in  the  reduction  steps  we  usually  eliminated  those  parameters  that  were  least 
significant;  i.o.,  those  chosen  least  often  in  the  subsets  of  the  previous  test; 
and,  generally,  the  lower  the  significance  of  a  parameter,  the  more  often  it  is 
missing  from  the  data  base.  It  is  concluded,  therefore,  that  the  decline  in  fore¬ 
cast  quality  in  Table  u  would  have  been  more  pronounced  had  all  parameters 
been  present  in  all  records.  This  proves  that  there  is  valuable  predictive  infor¬ 
mation  contained  in  at  least  some  of  the*  less  significant  parameters.  It  is 
emphasized  that,  perhaps  to  a  large  degree,  the  lower  significance  of  these 
parameters  is  due  only  to  their  frequent  absence  from  the  data  base. 

A  final  word  must  be  noted  regarding  the  combination  parameters.  Table  8 
indicates  that  a  number  of  these  new  parameters  have  been  selected  by  the  com¬ 
puter  program  as  significant  in  classifying  the  outcomes.  Due  to  the  complex 
intercorrelations  among  various  parameters,  however,  in  addition  to  possible 
variance  stabilization  effects  and  other  statistical  phenomena,  we  do  not  fully 
understand  the  true  significance  of  these  comoination  parameters,  (questions 
such  as  tins  probably  must  await  further  testing  on  data  bases  containing  fewer 
missing  parameters. 

l..'(  Texts  on  a  Kull\  Represented  Data  Base 

The  most  important  test  of  the  computer  forecast  is  achieved  in  the  case 
where  all  the  parameters  submitted  to  analysis  are  present  in  all  records  of  the 
data  base.  Such  a  test,  using  the  full  set  of  parameters,  is  impossible  with  the 
presently  available  data.  A  test  can  no  made  on  a  fully  represented  base,  how¬ 
ever,  if,  for  example,  only  eight  parameters  are  used,  and  we  are  willing  to 
accept,  a  reduced  base  of  3732  records,  of  which  only  22.32  remain  in  the  test  set. 
Such  a  test  (I)  was  performed,  and  the  results  are  shown  in  Tables  10,  11,  and  12. 

Test  I  shows  a  dramatic  improvement  in  the  MYDA/CL  computer  forecast 
in  all  scores,  while  the  .MYDA  and  comparison  forecasts  show  smaller  improve¬ 
ments.  These  improvements  occur  despite  the  somewhat  lower  flare  climatology 
that  applies  to  this  particular  test  set.  The  fact  that  the  comparison  (subjec¬ 
tive)  forecast  scores  are  higher  indicates  that  the  more  complete  observational 
coverage  during  this  sample  of  records  somehow  benefits  the  subjective  methods 
also. 

Due  to  the  reduced  number  of  records,  the  errors  associated  with  the  Test  I 
scores  are  about  50  percent  higher  than  those  stated  earlier.  Nevertheless,  there 
now  seems  no  question  that  the  MV  DA/ CL  forecast  is  superior  to  the  others. 
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Table  10.  Comparison  of  Forecasts  l  sing  a  Fully 
Represented  Data  Base- -Test  1  (1500-record  training  set) 


Largest  Event 

Largest  Event  Observed 

Tot  al 

Forecasted 

No  Flare 

C 

M  &  X 

Fo  recasts 

COMPARISON 

No  Flare 

17  54 

90 

8 

1852 

C 

193 

70 

24 

287 

M  &  X 

28 

31 

34 

93 

Total  Events 

1975 

191 

66 

2232 

MVDA 

No  Flare 

1707 

67 

10 

1784 

C 

232 

92 

24 

348 

M  &  X 

36 

32 

32 

100 

Total  Events 

1975 

191 

66 

2232 

MVDA/ CL 

No  Flare 

1829 

97 

14 

1940 

C 

145 

82 

33 

260 

M  &  X 

1 

12 

19 

32 

Total  Events 

197  5 

191 

66 

Table  11.  Parameters  Submitted  to  Analysis  and 
Their  Frequency  of  Selection  in  9  Subsets--Test  I 


Flares  Today 

9 

Mag.  Class 

8 

Spot  Class  2 

5 

Bright  Pts. 

9 

Mag.  Grad. 

8 

Spot  Class  1 

1 

Spot  Class  3 

8 

Spot  Dynam. 

8 
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I'ahle  12 


( 'on: pari. still  of  Forecast  Score.-.--  |  n.->t  1 


forecaster 

Event 

\ 

f 

A 

C  1  V.  <  iff  1 

<  iff  2 

COMPARISON 

No  flare 

94.7 

88. !! 

91.8 

88.  7 

: 

t: 

24.4 

:io.  r> 

'10.  7 

8.0  71.  7  84.  2  1  7.1 

1.0 

M  \  X 

-;g.  r; 

7 1 .  7 

4  4.  1 

o 

MY  1M 

No  Flare 

9  f> .  7 

KG.  4 

91.  1 

88.  7 

( ’ 

20.  4 

111.  2 

.‘17.  4 

8.  0  10.  2  82.  0  i  7.  “ 

2.1 

M  A.  \ 

;2.  •) 

7 

40.  2 

o 

M\  1 )  A  /  Cl. 

No  flare 

!,4 .  ", 

.-,2 

9  1.  7 

88.  '» 

C 

■ ;  1 .  •  ") 

■\2.  9 

.47.  2 

8.  0  .78.  4  80.  7  12.9 

0.  7 

M  \ 

:  ‘  .  i 

211.  d 

44.  1 

4.0 

:>.  conci.i  signs  \M) km:ommi.m>\iion> 

i'i)f  conclusions  of  thin  study  may  be  summarized  as  i'oilovrs: 

i.  The  standard  MY  DA  forecast  is  voi'v  similar  to  the  comparison  forecast 
used  in  this  stunt  in  terms  oi  overall  accuracy  and  bias  toward  overprediction. 

The  MYDA/CC  forecast  is  superior  overall  to  either  the  \1VI).\  or  the 
comparison  forecast,  and  is  biased  toward  underprediction. 

The  optimum  sire  for  file  training  set  is  probably  about  1.100  records 
for  the  climatoh  gies  that  prevailed  during  1977  and  1078. 

4.  f-'lares  l'oday  is  the  most  valuable  prediction  parameter  in  the  data 
base  used  here,  with  the  Bright  Points  parameter  a  verv  close  second. 

( >ther  important  parameters  are  Magnetic  Class,  Magnetic  Gradient,  Spot 
Class,  and  Sunspot.  Dynamics. 

7,  Combination  parameters,  although  their  role  is  not  fully  understood, 
seem  to  improve  forecast  scores. 

0.  Some  of  the  often  missing  parameters  (which  probably,  therefore,  only 
appear  to  be  less  significant  as  predictors)  contain  valuable  predictive  informa¬ 
tion.  Probable  candidates  include  Radio  Burst/ Sweep,  \eutral  Cine  Changes, 
Neutral  Cine  Complexity,  and  Emerging  Flux. 

The  MYDA/CC  procedure  may  be  capable  of  producing  forecasts  superior 
to  any  presently  available  using  conventional,  subjective  techniques.  It  has  been 
shown  that  its  skill  becomes  markedly  evident  when  complete  parameter  repre¬ 
sentation  is  achieved  in  the  data  lia.se.  On  the  basis  of  this,  we  predict  that  with 
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improvements  in  data  consistency,  as  well  as  the  inclusion  of  new,  objective 
parameters  in  the  future,  the  computer  forecast  scores  will  continue  to  improve. 


This  study  has  led  us  to  make  the  following  recommendations  concerning 
the  use  of  the  two  computer  forecasts: 

1.  Provide  a  flare  forecast  derived  from  MYDA/CL  for  those  customers 
who  cannot  tolerate  false  flare  alarms  (note  the  comparison  of  F  scores  in 
Table  12). 

2.  Provide  a  flare  forecast  derived  from  standard  MVDA  for  those  cus¬ 
tomers  who  need  to  be  forewarned  of  flares  as  often  as  possible  (compare  E 
scores  in  Table  12). 

3.  Improve  the  coverage  for  the  parameters  in  Table  1  that  are  deemed 
less  significant  by  virtue  of  their  frequent  absence  in  the  data  base. 

4.  Improve  the  objectivity  and  consistency  of  all  parameters. 


References 


Dixon,  'A  .  .1.  (ed. )  ( 1 9fiH >  Biomedical  Computer  Programs,  t  niversity  of 
California  Publications  in  Automatic  Computations,  No.  2,  l  niversity 
of  California  Press,  p.  214a. 

Seagraves,  P.  fl.  ( 1 97 2)  V ISC'  1JMD07 M  Stepwise  Discriminant  Analysis, 
Iniversity  of  British  Columbia  Computing  Centre  Documentation. 

Cooley,  W,  W.,  and  l.ohnes,  P.  H.  (1962)  Multivariate  Procedures  for  the 
Behavioral  Sciences,  Wiley,  New  York. 

Pachenbruch,  P.  \.,  and  Mickey,  M.  R,  (1968)  Technometrics,  10:1. 

\nclerson,  T.  i\ .  (195.8)  \n  Introduction  to  Multivariate  Statistical 
\nalysis,  Wiley, 

Kao,  C.  R.  (1974)  Advanced  Statistical  Methods  in  Biometric  Research, 
flafner. 

\  ecehia,  D.  I'.,  Caldwell,  (4.  A. ,  Tryon,  P.  \  . ,  and  .tones,  11.11. 

(1980)  in  Sol.  -Terrest.  Prod,  l’roc.  ,  Yol.  3,  11.  P.  Donnellv 
(eil.),  C-7fi. 

flirman,  J.U. ,  Neiciig,  ().  I'.  ,  Seagraves,  P.  IK,  1  Towers,  W  .  P . , 
and  Wiborg,  P.  H.  (1980)  in  Sol. -Terrest.  Pred.  Proc.,  Yol.  3, 

R.  F.  Donnelly  (ed.  ),  C-64. 

Simon,  P. ,  Smith,  J.B. ,  Ding,  Y.,  (Towers,  W  . ,  Guo,  Q. ,  Harvey, 

K.L.,  Hedeman,  H.,  Martin,  S.  F. ,  McKenna  Lawlor,  S.,  Pin,  Y., 
Neidig,  D. ,  Obridko,  Y.  Y. ,  Dodson  Prince,  IP,  Rust,  14.,  Speich, 

D. ,  Starr,  A.,  and  Stepanyan,  N.N.  (1980)  in  Sol. -Terrest.  Pred. 

Proc. ,  Yol.  2,  II.  F.  Donnelly  (ed.  ),  p.  287. 


